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1  Introduction 


The  development  of  residual  stress  is  strongly  influenced  by  processing  history.  Residual  stresses 
can  have  a  significant  effect  on  the  mechanical  performance  of  composite  structures  by  induc¬ 
ing  warpage  or  initiating  matrix  cracks  and  delaminations  [1,2, 3, 4].  Processing  concerns  associated 
with  thermosetting  composites  become  increasingly  important  for  components  of  appreciable  thick¬ 
ness.  The  most  common  problem  is  an  increase  in  temperature  resulting  from  the  resin  exothermic 
chemical  reaction  (polymerization)  that  may  raise  internal  temperatures  to  levels  inducing  mate¬ 
rial  degradation.  A  second  concern  is  the  complex  gradients  in  temperature  and  degree  of  cure 
accentuated  by  increased  thickness.  This  study  focuses  on  the  evolution  of  macroscopic  in-plane 
residual  stresses  in  thick  thermoset  laminates  resulting  from  complex  temperature  and  degree  of 
cure  gradients  that  develop  during  the  curing  process. 

Traditional  analyses  of  residual  stresses  in  thermosetting  composite  laminates  are  based  on 
thermal  expansion  mismatch  between  adjacent  plies,  a  uniform  temperature  difference  between 
the  cure  temperature  and  ambient  and  no  stress  development  prior  to  completion  of  the  curing 
process  [3, 4, 5,6].  These  analyses  have  been  quite  successful  in  predicting  residual  stresses  in  thin 
laminates,  where  a  uniform  through-the-thickness  temperature  distribution  assumption  is  justified. 
Such  an  approach,  however,  is  not  appropriate  for  predicting  the  process-induced  stresses  in  thick 
thermoset  composite  laminates  where  complex  temperature  and  degree  of  cure  gradients  develop 
during  the  curing  process  [7,8,9,10,11,12,13].  As  we  demonstrate,  these  temperature  and  degree  of 
cure  gradients  represent  important  mechanisms  that  contribute  to  the  development  of  stress  and 
deformation  in  thick-section  composite  laminates. 

These  mechanisms  are  quite  similar  to  those  which  govern  the  development  of  residual  stresses 
encountered  in  the  manufacture  of  tempered  glass  [14,15,16,17].  Stresses  result  from  the  interactions 
between  spatially  varying  thermal  contractions  and  viscoelastic  material  response  induced  by  severe 
temperature  gradients  which  develop  during  the  quenching  process.  Early  treatment  of  this  type  of 
stress  development  due  to  viscoelastic  material  response  in  the  presence  of  thermal  gradients  were 
based  on  the  time- temperature  superposition  principle  of  Theologically  simple  material  behavior  [18]. 
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This  approach  has  been  applied  to  study  the  macroscopic  in-plane  residual  stress  development  under 
the  rapid  cooling  process  of  an  isotropic  epoxy  plate  [19]  and  more  recently  to  the  quenching  of 
thermoplastic  matrix  composite  laminates  [20]. 

Complex  temperature  and  degree  of  cure  gradients  which  develop  in  thick-section  thermoset¬ 
ting  laminates  during  the  curing  process  induce  spatially  varying  material  response  similar  to  that 
which  induces  viscoelastic  stress  development  in  rapid  cooling  and  quenching  processes.  Bire¬ 
fringence  patterns  have  shown  the  effect  of  temperature  and  degree  of  cure  gradients  on  the  de¬ 
velopment  of  residual  stresses  in  thick  thermoset  castings  [7],  In  a  series  of  papers,  Levitsky  and 
Shaffer  [21,22,23]  have  studied  the  development  of  residual  stresses  including  temperature  gradients 
and  spatially  varying  chemical  hardening  effects  on  mechanical  properties  in  isotropic  thermosetting 
materials  accentuated  by  the  exothermic  chemical  reaction.  These  studies  indicate  that  tempera¬ 
ture  and  degree  of  cure  gradients  significantly  influence  stress  development  prior  to  a  fully  cured 
state  and  consequently  should  be  accounted  for  in  studying  the  evolution  of  process-induced  stress 
and  deformation  in  thick-section  thermosetting  composite  laminates. 

Another  important  mechanism  contributing  to  stress  development  not  considered  in  the  analyses 
cited  [21,22,23]  is  the  volumetric  shrinkage  of  the  thermosetting  resin  associated  with  the  cross¬ 
link  polymerization  reaction.  The  importance  of  chemical  shrinkage,  which  is  an  additional  source 
of  internal  loading,  can  be  appreciated  by  considering  the  shrinkage  occurring  within  the  interior 
region  of  a  laminate  that  is  constrained  by  a  fully  cured  exterior  region.  A  similar  situation  occurs 
on  the  micromechanics  level  where  resin  shrinkage  during  cure  is  constrained  by  the  reinforcement 
phase  [24]. 

As  a  thermosetting  resin  cures  its  material  characteristics  change  dramatically,  transitioning 
from  the  behavior  of  a  viscous  liquid  (low  stiffness),  in  its  uncured  state,  to  a  viscoelastic  or 
elastic  solid  (high  stiffness),  in  its  fully  cured  state.  Significant  reduction  in  specific  volume  (cure 
shrinkage)  associated  with  the  cross-link  polymerization  reaction  also  occurs  during  cure.  These 
changes  in  the  resin  directly  influence  the  effective  mechanical  properties  and  cure  shrinkage  strains 
of  the  composite  during  cure. 

In  the  present  investigation,  the  curing  process  of  the  resin  is  separated  into  three  distinct 
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regions  for  descriptive  purposes  as  shown  in  Figure  1.  In  Region  I,  the  resin  is  in  the  B-stage 
condition,  fully  uncured  and  assumed  to  behave  as  a  viscous  fluid  (negligible  stiffness).  Region  II 
denotes  the  curing  stage  of  the  resin,  where  a  significant  increase  in  stiffness  (chemical  hardening) 
and  a  reduction  in  specific  volume  (chemical  shrinkage)  occurs.  In  this  region  the  mechanical 
properties  of  the  resin  are  governed  inherently  by  competing  mechanisms  between  chemical-kinetic 
hardening  and  viscoelastic  relaxation  phenomena.  Region  III  marks  the  end  of  the  curing  process 
and  no  further  chemical  shrinkage  occurs.  In  Region  III  the  resin  exhibits  traditional  viscoelastic 
behavior  at  elevated  temperatures  and  approaches  elastic  behavior  at  lower  temperatures.  Thermal 
expansion  is  the  only  mechanism  contributing  to  changes  in  specific  volume  in  Region  III. 

This  work  focuses  on  the  influence  the  curing  process  has  on  the  macroscopic  in-plane  stress 
and  deformation  development  in  thick  thermosetting  composite  laminates.  A  one-dimensional  cure 
simulation  model,  based  on  an  incremental  transient  finite  difference  formulation  that  accounts 
for  thermal  and  chemical  interactions,  is  developed  to  predict  complex  temperature  and  degree  of 
cure  gradients  that  develop  during  cure.  Material  models  are  proposed  to  describe  the  modulus  and 
shrinkage  of  the  resin  during  cure.  Mechanical  properties  of  the  fiber  phase  are  assumed  independent 
of  cure.  A  micromechanics  model,  employing  the  resin  and  fiber  constituent  behavior,  is  used  to 
evaluate  the  instantaneous  spatially  varying  mechanical  properties,  thermal  expansion  and  chemical 
shrinkage  strains  within  the  composite  laminate  as  a  function  of  temperature  and  degree  of  cure. 
Process-induced  stress  and  deformation  predictions  are  based  on  an  incremental  laminated  plate 
theory  model  that  includes  temperature  gradients,  spatially  varying  cure  dependent  mechanical 
properties,  thermal  expansion  and  chemical  shrinkage  strains.  The  incremental  laminated  plate 
theory  model  and  cure  simulation  analysis  are  coupled  enabling  the  prediction  of  macroscopic 
in-plane  process-induced  stress  and  deformation  in  thick  composite  laminates  during  the  curing 
process  to  be  made. 

The  methodology  developed  for  predicting  process-induced  stress  and  deformation  in  thick 
sections  is  coded  into  a  FORTRAN  computer  program.  The  program  is  employed  in  parametric 
studies  to  quantify  the  influence  of  processing  on  residual  stress  development  in  thick  glass/polyester 
and  graphite/epoxy  thermoset  laminates.  Model  predictions  of  cure  dependent  epoxy  modulus 
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Figure  1:  Resin  modulus  and  chemical  shrinkage  during 


cure 
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and  curvature  in  unsymmetric  graphite/epoxy  laminates  are  correlated  with  experimental  data. 
Processing  history  (autoclave  temperature  cycle),  laminate  thickness  and  resin  cure  shrinkage  are 
shown  to  have  a  profound  influence  on  the  residual  stress  development  in  thick  sections.  The 
influence  of  unsymmetric  curing  (curing  the  laminate  from  one  side)  and  laminate  stacking  sequence 
on  process-induced  stress  development  is  also  investigated.  Results  are  presented  which  demonstrate 
that  the  traditional  assumption  of  a  stress-free  laminate  temperature  is  not  generally  appropriate 
for  predicting  residual  stresses  in  thick  themoset  composite  laminates. 

2  Analysis 

2.1  Cure  Simulation 

The  one-dimensional  cure  simulation  analysis  presented  herein  was  developed  as  a  simplification 
of  the  authors’  full  two-dimensional  anisotropic  cure  simulation  analysis  for  complex  shaped  cross- 
sectional  geometries  [9,13].  A  one-dimensional  analysis  enables  isolation  of  through-the-thickness 
processing  effects  on  a  fundamental  level  without  the  extra  computational  effort  required  in  a  two- 
dimensional  analysis.  The  solution  strategy  utilized  in  the  cure  simulation  analysis  is  consistent 
with  the  treatment  made  by  previous  investigators  [25,26]-  The  principle  governing  equation  utilized 
in  the  analysis  is  the  one-dimensional  Fourier’s  Heat  Conduction  Equation: 

fl2rp  dT 

^Jrkz'dz^~pCp~di  for  T(z>t)  in  (°  <  z  <  l)  (1) 

The  internal  heat  generation  term,  g,  accounts  for  the  exothermic  chemical  reaction  characteristic 
of  thermosetting  systems.  The  parameters  kz,  p  and  cp  are  the  effective  interlaminar  (through- 
the-thickness)  thermal  conductivity,  density  and  specific  heat  of  the  composite,  respectively.  These 
thermal  properties  are  assumed  constant  during  the  curing  process  [7].  Total  laminate  thickness  is 
t.  Absolute  temperature  and  time  are  denoted  T  and  t ,  respectively.  The  z—  direction  is  normal  to 
the  in-plane  dimension  of  a  laminate.  Thermal  properties  for  glass/polyester  and  graphite/epoxy 
used  in  this  study  are  summarized  in  Table  1  [25,26,27]. 

A  generalized  temperature  boundary  condition  formulation  was  employed  allowing  for  either 
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P  [kg/™3] 

cp  [kJ/(W  °C)] 

kz  [kW/(m  •  °C)] 

Glass/Polyester 

Graphite/Epoxy 

1.89xlOa 

1.52xl03 

1.26 

9.42xl0-1 

2.163xl0~4 

4.457x10““ 

Table  1:  Thermal  properties  for  glass/polyester  and  graphite/epoxy. 


convective,  insulated  or  prescribed  temperature  boundary  conditions  on  the  laminate  surface: 

+  h'fj{T,  -  T(t))  =  0 

at  z  =  0  and  z  =  £  (2) 

The  temperature  and  normal  derivative  of  temperature  on  the  laminate  surface  are  denoted  T,  and 
dTgfdz ,  respectively.  The  coefficients  keff  and  hejj  represent  the  effective  thermal  conductivity  and 
convective  heat  transfer  coefficient  on  the  laminate  surface,  respectively.  The  autoclave  temperature 
cure  cycle  is  incorporated  into  the  cure  simulation  through  T(t).  Prescribed  temperature  boundary 
conditions  (keff= 0  and  hejj=l)  were  employed  in  all  simulations  presented  in  this  work  to  eliminate 
the  added  complexity  of  interpreting  the  influence  of  convection  on  the  results  [9,13]. 

The  internal  heat  generation  term  in  equation  (1),  q,  represents  the  instantaneous  heat  liberated 
per  unit  volume  of  material  from  the  cross-link  polymerization  reaction: 

da 

9  =  PHt  (3) 

The  heat  of  reaction,  Hr,  is  the  total  heat  liberated  for  complete  cure  and  da/dt  is  the  instantaneous 
cure  rate.  The  degree  of  cure  at  any  time  is  defined  in  terms  of  the  instantaneous  cure  rate  through 
an  integral  representation: 

Jft  da 

[  nJI  <4> 

The  complete  description  of  the  cure  kinetics  for  the  composite  includes  the  total  heat  of  reac¬ 
tion  and  a  description  of  the  rate  of  reaction  as  a  function  of  temperature  and  degree  of  cure.  The 
instantaneous  reaction  rate  is  required  to  compute  the  heat  generation  (equation  (3))  and  degree  of 
cure  (equation  (4))  during  the  curing  process.  Both  the  total  heat  of  reaction  and  the  reaction  rate 
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expression  are  typically  characterized  empirically  with  isothermal  Differential  Scanning  Calorime¬ 
try  (DSC)  techniques  [27].  Reaction  rate  expressions  for  glass/polyester  and  graphite/epoxy  are 
different  in  form  due  to  the  inherent  differences  in  the  overall  order  of  the  reaction  kinetics. 

The  glass/polyester  composite  contains  CYCOM  4102  polyester  resin,  manufactured  by  the 
American  Cyanamid  Corporation,  and  is  reinforced  with  E-glass  fibers  (54%  by  volume).  The 
reaction  rate  expression  for  glass/polyester  is  second-order  overall  in  the  sense  that  (mc  +  nc  = 
2),  [27]: 

~  =  Acexp(-AEc/RT)am‘(l  -  a)n‘  (5) 

The  parameter  R  denotes  the  universal  gas  constant.  The  exponents  mc  and  nc,  the  pre-exponential 
coefficient,  Ac,  the  activation  energy,  A Ec,  and  the  total  heat  of  reaction,  Hr>  for  glass/polyester 
are  listed  in  Table  2. 

The  graphite/epoxy  composite  contains  Hercules  Corporation’s  3501-6  epoxy  resin  and  is  rein¬ 
forced  with  unidirectional  AS4  graphite  fibers  (67%  by  volume).  The  reaction  rate  expression  for 
graphite/epoxy  follows  a  markedly  different  form  [25,26]: 


=  (*1  +  *2<>)(1  ~  «)(0.47  -  a) 

for 

(a  <  0.3) 

da 

—  =  t3(i  -  <*> 

for 

(a  >  0.3) 

(6) 

The  parameters  ifcx,Jb2  and  £3  are  defined  by  the  Arrhenius  rate  expressions 

k\  =  A\exp{— AE\  /  RT) 

k2  =  A2exp(-AE2/RT)  (7) 

£3  =  A3exp(—AE3/ RT) 

The  pre-exponential  coefficients  A\,  A2  and  A3,  the  activation  energies,  A E\,  AE2  and  AE3,  and 
the  total  heat  of  reaction  for  graphite/epoxy  are  summarized  in  Table  2  [25,26,27]. 

Initial  conditions  for  temperature  and  degree  of  cure  are  mathematically  represented  by 

T(z)  =  Ti  for  (0  <  z  <  £)  at  t  =  0 

a(z)  =  a i  for  (0  <  z  <  t)  at  t  =  0  (8) 
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Glass/ Polyester 

mc 

nc 

Ac[min.-1] 
A£c[J/mol] 
Hr  [k  J  /kg] 

0.524 

1.476 

3.7xl022 

1.674xl05 

77.5 

Graphite/Epoxy 

Aj[min.-1] 

Ajfmin.-1] 

A3[min.-1] 

A  Ei  [J/mol] 
A  £2  [J /mol] 
A  £3  [J/mol] 
tfr[kJ/kg] 

2.102xl0a 

-2. 014x10* 

1.960xl05 

8.07x10“ 

7.78x10“ 

5.66x10“ 

198.9 

Table  2:  Cure  kinetic  parameters  for  glass/polyester  and  graphite/epoxy. 


T;  and  a;  are  the  prescribed  uniform  initial  temperature  and  degree  of  cure  distributions  within 
the  laminate  taken  to  be  the  ambient  temperature  and  zero,  respectively. 

An  incremental  transient  finite  difference  technique  is  employed  to  solve  the  governing  equations, 
boundary  and  initial  conditions  that  define  the  cure  simulation  in  question.  Transient  tempera¬ 
ture  and  degree  of  cure  distributions  (through-the-thickness)  are  predicted  as  a  function  of  the 
thermal  properties,  chemical-kinetic  parameters,  initial  and  boundary  conditions  including  the  au¬ 
toclave  temperature  cure  cycle.  The  laminate  is  discretized  in  the  thickness  dimension  (2)  with 
a  one-dimensional  finite  difference  nodal  grid.  The  spatial  and  time  derivatives  in  the  governing 
temperature  equation  (1)  and  the  boundary  condition  (2)  are  replaced  by  their  appropriate  finite 
difference  approximations  [9,13].  Transient  temperature  distributions  are  solved  from  the  result¬ 
ing  set  of  finite  difference  equations  using  the  Alternating  Direction  Explicit  (ADE)  Method  [28]. 
Degree  of  cure  distributions  are  computed  from  a  finite  sum  approximation  of  the  integral  rep¬ 
resentation  given  by  equation  (4).  The  solution  strategy  of  the  cure  simulation  marches  out  the 
autoclave  temperature  cure  cycle  until  its  completion.  In  this  manner,  the  transient  temperature 
and  degree  of  cure  distributions  within  the  laminate  are  generated  as  a  function  of  the  processing 
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history.  Details  of  the  solution  strategy  are  documented  elsewhere  [9,13]. 


2.2  Material  Models 

Two  material  models  are  proposed  to  describe  the  behavior  of  the  thermoset  resin  during  cure. 
These  models  describe  the  mechanical  property  changes  and  volumetric  cure  (chemical)  shrinkage 
associated  with  the  cross-linked  polymerization  reaction.  Property  changes  and  cure  shrinkage  of 
the  resin  during  cure  directly  influence  the  homogeneous  mechanical  properties  and  process-induced 
strains  in  the  composite.  This  influence  is  quantified  with  a  micromechanics  mode!  for  unidirectional 
continuous  fiber  reinforced  composites.  Mechanical  property  and  process-induced  strain  gradients 
through  the  laminate  thickness  represent  important  mechanisms  contributing  to  macroscopic  in¬ 
plane  stress  and  deformation  development  during  cure.  The  material  models  that  describe  the  resin 
behavior  during  cure,  coupling  the  cure  simulation  and  stress  analysis,  are  discussed  below. 

2.2.1  Cure  Dependent  Resin  Modulus 

The  resin  modulus  model  describes  the  mechanical  properties  of  the  resin  during  cure.  The  resin 
modulus  is  strongly  cure  dependent,  influenced  by  the  kinetic-viscoelastic  interactions  successfully 
modeled  by  Dillman  and  Seferis  [29].  While  their  model  was  rigorous,  independent  evaluation  of 
the  kinetic  and  viscoelastic  parameters  required  extensive  data  reduction  procedures.  In  addition, 
model  predictions  outside  the  temperature  range  cover  by  Dillman  and  Seferis  were  not  reliable. 
Consequently,  a  more  convenient  a-mixing  rule  model  [30]  is  proposed  here  to  describe  the  kinetic- 
viscoelastic  behavior  of  the  resin  modulus  during  cure.  The  instantaneous  isotropic  resin  modulus, 
denoted  Em,  is  expressed  explicitly  in  term  of  degree  of  cure: 

Em  =  (1  -  a)E°m  +  +  7a(l  -  a)(E%  -  E°m)  (9) 

where  (— 1  <  7  <  1) 

The  parameters  and  E™  are  the  assumed  fully  uncured  (Region  I)  and  fully  cured  (Region 
III)  temperature  dependent  resin  moduli,  respectively.  The  term  7  is  introduced  to  quantify  the 
competing  mechanisms  between  stress  relaxation  and  chemical  hardening  [29].  Increasing  7  physi- 
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Table  3:  Resin  characteristics  during  cure. 


cally  corresponds  to  a  more  rapid  increase  in  modulus  at  lower  degree  of  cure  before  asymptotically 
approaching  the  fully  cured  modulus.  Note  that  this  model  is  extremely  sensitive  to  the  kinetics 
of  the  resin  and,  therefore,  the  non-isothermal  process  history  (see  equations  (5)  and  (6)).  Results 
presented  in  this  study  assume  E and  £“  are  constant  and  7  =  0. 

The  instantaneous  resin  shear  modulus  during  cure  is  based  on  the  isotropic  material  relation: 

C  — 

—  2(1  +  vm) 

Poisson’s  ratio  of  the  resin,  um ,  is  assumed  constant  during  cure.  As  this  investigation  focuses  on 
the  influence  that  the  curing  process  has  on  stress  and  deformation  development,  all  mechanical 
properties  of  the  resin  are  assumed  constant  once  cure  is  complete.  Therefore,  the  analysis  predic¬ 
tions  for  this  assumption  may  provide  upper  bounds  on  process-induced  stresses,  since  no  stress 
relaxation  is  considered. 

The  fully  uncured  and  fully  cure  resin  moduli  for  glass/polyester  and  graphite/epoxy  used  in 
this  study  are  listed  in  Table  3  [31].  Uncured  moduli  (£'£l)  are  chosen  arbitrarily  small  (negligible 
stiffness)  while  the  fully  cured  moduli  ( E %?)  are  representative  of  typical  room  temperature  values. 

In  Figure  2,  the  mode!  of  Levitsky  and  Shaffer  [21,22,23]  is  correlated  with  the  present  formula¬ 
tion  for  the  chemical  hardening  of  a  polyester  resin  (see  Tables  2  and  3)  during  isothermal  cure  at 
100  °C.  The  present  model  assumes  the  Poisson  ratio  is  constant  while  the  latter  model  assumes  the 
bulk  modulus  is  constant.  Both  models  predict  nearly  identical  results  for  modulus  development 
during  cure.  Deviation  in  Poisson  ratio  is  noted  at  low  degree  of  cure  while  the  models  converge 
rapidly  as  the  resin  cures  (see  Figure  2).  The  macroscopic  composite  properties,  process-induced 
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MODULUS  (MPa) 


TIME  (MINUTES) 


Figure  2:  Influence  of  the  constant  poisson  ratio  assumption  on  modulus  development  during  cure. 
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POISSON  RATIO 


strains  and  residua]  stresses  are  not  significantly  influenced  by  this  difference  in  the  Poisson  ratio 
of  the  resin  during  cure. 

2.2.2  Cure  Dependent  Resin  Chemical  Shrinkage 

A  second  material  model  is  proposed  to  describe  the  volumetric  chemical  shrinkage  of  the  resin 
during  cure.  Resin  shrinkage  only  occurs  during  the  curing  process  and  ceases  once  cure  is  complete. 
Chemical  resin  shrinkage  induces  significant  macroscopic  strains  in  the  composite,  representing  an 
important  source  of  internal  loading  in  thick-section  laminates  in  addition  to  the  traditionally 
recognized  thermal  expansion  strains. 

The  volumetric  change  of  a  cubic  volume  element  of  dimension  / j  by  /2  by  I3  can  be  expressed 
in  terms  of  its  overall  dimensions  and  the  finite  dimensional  changes  in  three  principle  directions, 
A/i,  A/2,  A/3,  as 

AV'  =  /1A/2/3  +  A/1/2/3  +  A/1A/2/3  +  /1/2A/3  + 

/1A/2A/3  4-  A/1/2AI3  -f  A/1A/2A/3  (11) 

An  associated  change  in  specific  volume,  Ad,  can  be  defined  in  terms  of  the  principle  strain  com¬ 
ponents, 

AV 

Ad  =  =  €j  +  «2  -f  £3  +  fi«2  +  C1C3  +  £2f3  +  «i£2f3  (12) 

Assuming  a  uniform  strain  contraction  for  all  principle  strain  components,  the  incremental 
isotropic  shrinkage  strain,  Aer,  of  a  unit  volume  element  of  resin  resulting  from  an  incremental 
specific  volume  resin  shrinkage,  Avr,  becomes1 

Acr  =  ^1  +  Ad,  -  1  (13) 

The  incremental  volume  resin  shrinkage  is  based  on  an  incremental  change  in  degree  of  cure,  Aa, 
and  the  total  specific  volume  shrinkage  of  the  completely  cured  resin,  vjhl  through  the  following 
expression: 

Avr  =  Aa-vJh  (14) 

’Exact  expression  suggested  by  Joseph  M.  Santiago,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground,  MD. 
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Typical  total  volumetric  resin  cure  shrinkage  values  for  glass/polyester  and  graphite/epoxy  are 
listed  in  Table  3. 

2.2.3  Cure  Dependent  Composite  Mechanical  Properties 

The  effective  homogeneous  unidirectional  mechanical  properties  of  each  individual  lamina  within 
the  composite  laminate  are  computed  each  time  increment  during  the  cure  simulation.  Lamina 
properties  are  highly  dependent  on  the  fiber  and  resin  constituent  properties,  and  fiber  volume 
fraction.  The  mechanical  properties  of  the  resin  (except  poisson  ratio  and  thermal  expansion 
which  are  assumed  constant)  vary  according  to  the  material  models  presented  above.  Mechanical 
properties  of  the  fiber  are  assumed  constant  and  independent  of  cure.  During  cure,  significant 
composite  property  gradients  develop  as  a  result  of  the  temperature  and  degree  of  cure  gradients 
associated  with  thick  laminate  processing. 

The  doubly-embedded  self-consistent  field  micromechanics  model  [32]  is  used  to  compute  the 
instantaneous  transversely  isotropic  mechanical  properties  and  thermal  expansion  coefficients  of 
the  lamina.  The  micromechanics  equations  employed  in  this  investigation  are  documented  in 
the  Appendix.  The  fiber  and  resin  constituent  mechanical  properties  for  glass/polyester  and 
graphite/epoxy  are  summarized  in  Table  4  [31].  The  1-direction  is  coincident  with  the  direction  of 
fiber  reinforcement.  The  2-direction  (in-plane)  and  3-direction  (out-of-plane)  are  both  perpendic¬ 
ular  to  the  1-direction.  The  graphite  fiber  is  assumed  transversely  isotropic  while  the  glass  fiber  is 
assumed  isotropic.  Note  the  strong  resin  property  dependency  on  degree  of  cure  and  thus  on  the 
processing  history. 

2.2.4  Incremental  Composite  Chemical  Shrinkage  Strain 

Effective  chemical  shrinkage  strains  in  the  composite  are  also  computed  according  to  the  microme¬ 
chanics  model  and  are  based  on  the  fiber  and  resin  mechanical  properties,  chemical  resin  shrinkage 
strain  and  fiber  volume  fraction.  The  in-plane  principle  lamina  1-direction  (longitudinal)  and  2- 
direction  (transverse)  chemical  shrinkage  strain  increments,  denoted  &e\h  and  Ae^,  respectively, 
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Property 

Glass 

Polyester 

Graphite 

Epoxy 

Ei  [MPa] 

7.308x10“ 

eqn.(  10) 

2.068xl05 

eqn.(10) 

£2  [MPa] 

7.308x10“ 

eqn.(10) 

2.068x10“ 

eqn.(10) 

"12 

0.22 

0.40 

0.20 

0.35 

"13 

0.22 

0.40 

0.20 

0.35 

"23 

0.22 

0.40 

0.50 

0.35 

Gn  [MPa] 

2.992xl04 

eqn.(10) 

2.758x10“ 

eqn.(10) 

G13  [MPa] 

2.992x10“ 

eqn.(10) 

2.758x10“ 

eqn.(10) 

G23  [MPa] 

2.992xl04 

eqn.(10) 

6.894xl03 

eqn.(10) 

a,  [1/‘C] 

5.04xl0-6 

7.20xl0~5 

-9.00xl0-7 

5  76xl0-5 

a2  [1/°C] 

5.04xl0-6 

7.20xl0~5 

7.20xl0-6 

5.76xl0~5 

Table  4:  Fiber  and  resin  constituent  mechanical  properties. 


are  computed  over  each  time  increment  in  the  cure  simulation.  Specifically,  these  strains  are  based 
on  the  expansional  strain  relations  of  the  micromechanics  model  utilizing  the  instantaneous  me¬ 
chanical  properties  of  the  fiber  and  resin,  the  resin  chemical  shrinkage  strain  increment,  a  zero 
stress-free  shrinkage  strain  in  the  fiber  and  the  fiber  volume  fraction. 

2.2.5  Incremental  Composite  Thermal  Expansion  Strain 

Incremental  thermal  expansion  strains  are  also  computed  over  each  time  increment  during  the 
cure  simulation.  They  are  based  on  the  lamina  temperature  increment  between  two  consecutive 
time  steps,  AT,  and  the  instantaneous  effective  longitudinal  and  transverse  thermal  expansion 
coefficients,  aj  and  02,  respectively.  The  incremental  longitudinal  and  transverse  strain  increments 
in  each  lamina  are  defined  in  the  usual  manner: 

=  a-]  ■  AT 

A4*=o2-A  T  (15) 
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2.2.6  Total  Incremental  Composite  Process-Induced  Strain 


The  total  stress-free  macroscopic  process-induced  lamina  strain  increment  (over  a  single  time  step) 
is  computed  by  superimposing  the  thermal  and  chemical  strain  contributions.  The  total  process- 
induced  lamina  strain  increments  in  the  longitudinal  and  transverse  directions  become 

AeJ  =  Ae\h  +  Acf* 

=  Ae?  +  Aef  (16) 

Both  chemical  (shrinkage)  and  thermal  (expansional)  strains  are  inherently  elongational,  con¬ 
sequently  no  shear  strains  in  the  principle  material  coordinate  system  are  required  in  the  analysis. 
Note,  however,  that  process-induced  shear  strains  will  develop  in  the  global  laminate  system  when 
other  than  cross-ply  [0/90]  laminate  constructions  are  considered. 

2.3  Process-Induced  Stress  Modeling 

Classical  laminated  plate  theory  forms  the  theoretical  basis  for  the  stress  analysis  presented  herein  [33, 
34].  An  incremental  Hooke’s  law  formulation  is  employed  which  utilizes  the  instantaneous  effective 
composite  properties  and  incremental  process- induced  lamina  strain  increments  as  input.  Vari¬ 
ous  incremental  laminated  plate  theory  solution  strategies  and  convergence  criteria  are  discussed 
elsewhere  [35].  The  theory  and  methodology  employed  in  this  study  are  discussed  below. 

Prior  to  the  stress  and  deformation  calculations  for  a  single  time  step  increment,  the  cure  sim¬ 
ulation  is  performed  to  yield  the  temperature  and  degree  of  cure  distributions  within  the  laminate. 
From  these  distributions,  the  instantaneous  effective  mechanical  properties  and  incremental  process- 
induced  strains  in  the  laminate  are  computed  based  on  the  material  models  discussed  above.  Stress 
and  strain  increments  are  then  computed  from  the  plate  theory  analysis  over  each  time  step  in  the 
cure  simulation.  This  procedure  is  repeated  for  each  time  step  increment  of  the  cure  simulation.  By 
maintaining  a  cumulative  sum  of  the  stress  and  strain  increments,  a  complete  description  of  their 
transient  development  during  the  curing  process  and  subsequent  cool-down  to  uniform  ambient 
temperature  is  obtained. 
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The  laminate  is  discretized  with  nodal  points  in  the  2-direction,  consistent  with  the  finite 
difference  grid  employed  in  the  cure  simulation  analysis.  Mesh  refinement  is  established  such  that 
a  ply  is  defined  between  each  pair  of  nodes  in  the  finite  difference  grid,  uniformly  spaced  a  distance 
A z.  The  temperature  and  degree  of  cure  of  a  ply  is  taken  as  the  arithmetic  average  of  the  cure 
simulation  predicted  values  at  the  top  and  bottom  nodes  of  that  ply. 

A  fixed  global  laminate  coordinate  system  is  established  where  x-  and  y-  are  the  in-plane  direc¬ 
tions  and  z-  is  the  out-of-plane  direction.  The  stress  analysis  presented  here  addresses  macroscopic 
in-plane  stress  development  induced  by  processing  gradients  in  the  out-of-plane  direction.  The 
principle  (1,  2)  coordinate  system  of  each  ply  can  be  orientated  any  angle  6  with  respect  to  the 
fixed  (x,  y)  global  coordinate  system  of  the  laminate.  Therefore,  the  analysis  presented  herein 
accommodates  laminates  of  arbitrary  stacking  sequence. 

The  effective  in-plane  force  and  moment  resultants  are  first  calculated  for  the  time  step  incre¬ 
ment.  The  incremental  in-plane  force  resultants,  ANj ,  ANj  and  ANjy,  and  incremental  moment 
resultants,  AMj ,  AMy  and  AMjy,  are  given  by: 

(AATJ,  AA/J)  =  +  QuAel  +  Q\  A^XM) 

{ANj ,  AMy  )  =  ELi +  Q\^Tv  +  Q£6AeJy)(l,2)  (17) 

(AN^y,AMjy)  =  Z?=1(Q}6AJx  +  Q\  AJy  +  Q^A^Kl,*) 

The  total  number  of  plies  in  the  laminate  is  N.  The  matrix,  Qjj,  is  the  transformed  plane  stress 
stiffness  matrix  of  the  kth  ply  defined  in  the  global  ( x,y )  laminate  coordinate  system  in  the  tradi¬ 
tional  manner  [33,34]: 

Qn  =  ™4Qn  +  2m2n2(Qi2  4-  2^66)  +  n4Q22 
Qu  =  m2n2(QH -(- Q22  -  4Qe6)  +  (”i4  +  n4)Qi2 

Q22  =  n4Qn  +  2m2n2(Qi2  -f  2Qee)  4-  m4Q22 

Qi6  =  *ri3n((5ii  -  Q\2  ~  2^66)  +  mn3(Q\2  —  Q22  +  2^66) 

Q26  =  mn3(Qu  —  Q12  -  2Qee)  +  rn3n{Qn  -  Q22  +  2Qee) 

Q&&  =  ™2n2(Qu  +  Q22  -  2Qi2  -  2<566)  +  (m4  +  n4)<566  (18) 
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where  m  =  cosd  and  n  =  sind  and  the  nonzero  stiffness  constants  in  the  principle  material  coordi¬ 
nate  system  are  given  by 

<2..  = 

(1  -  ^12^21) 

V12  =  7: - r 

(1  —  1^12^21) 

Q27  =  7-j - - - 7 

(1  -  1^1 2*^21  ) 

C?66  =  G\2  (19) 


The  stress-free  incremental  ply  strains,  Aex ,  Aey  and  AcJy,  defined  in  the  global  laminate  coor¬ 
dinate  system,  are  based  on  the  incremental  stress-free  process-induced  ply  strains  in  the  principle 
lamina  coordinate  system  (equation  16)  through  the  second-order  tensoral  transformation 


m  n  mn 

„2  _ 


n*  m"  — mn 

—  2  mn  2mn  (m2  —  n2) 


The  effective  incremental  laminate  in-plane  strains,  Ae°,  Ae°  and  Ae°  ,  and  curvatures,  Akx, 


Aks  and  A/cXy  are  then  evaluated: 


A  Nj 

ANl 

am/ 

A  Ml 

A  AA* 


Ae®  an  <*12  <*16  612  b\G  ANX 

Ae°  a]  2  022  a26  £>12  622  £>26  AAf/ 

(  Ae°y  _  ai6  a26  616  £>26  £>66  AiVxy  ^  , 

Akx  £>n  £>12  £>16  dn  d di6  AMT  '  1 

A/Cv  612  £>22  £>26  d\2  ^22  <^26  AM/ 

.  AfCXy  L  £>16  £>26  £>66  ^16  d2 6  ^66  J  , 

The  matrices  [a,y],  [6jj]  and  [djj]  are  the  instantaneous  effective  laminate  compliance  coefficients 
computed  in  the  traditional  manner  [33,34]  with  spatial  variation  of  instantaneous  mechanical  ply 
properties  taken  into  account. 

Once  the  incremental  in-plane  strains  and  curvatures  are  determined,  the  laminate  ply  strain 
increments,  Aex,  Aev  and  Aexv  are  computed  through  the  classic  strain-displacement  relations  [33, 
34]: 

Aex  =  Ae°  +  ZkAnx 

Acy  =  Ac"  +  Zk  A«y  (22) 


Afxy  =  Ac"  +  zk  Akx 
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Here  Zk  is  the  distance  from  the  laminate  mid-plane  to  the  center  of  the  kth  ply. 

The  incremental  ply  stresses  are  based  on  the  difference  between  the  laminate  ply  strain  incre¬ 
ments  and  the  stress-free  process-induced  ply  strain  increments  through  the  following  expression: 

Acr*  =  Q\(Aex-Ael)k  +  Q\(Aey-Ael)k  +  Q]6(&exy-Aely)k 
A  crj  =  q}2(  Aex  -  AcTx)k  +  Aey  -  AeTy)k  +  Qk6(  Aexy  -  AeTxy)k  (23) 

A  aky  =  Qk6(Aex-Ael)k  +  Qke(Aey-Ael)k  +  Qk6(Aexy-Aely)k 

Transient  process-induced  stress  and  strain  distributions  through  the  laminate  thickness  are  ob¬ 
tained  from  the  cumulative  sum  of  the  increments  computed  at  each  time  6tep  during  the  cure 
simulation.  The  governing  equations  for  the  cure  simulation,  material  models  and  stress  analysis 
were  coded  into  a  FORTRAN  computer  program.  A  flow  diagram  summarizing  the  key  steps  in 
the  solution  strategy  is  presented  in  Figure  3.  As  an  example  of  the  computational  efficiency  of  the 
computer  simulation,  a  2.54  cm  thick  laminate  cured  for  six  hours,  consisting  of  30  nodes  and  10000 
time  steps,  takes  approximately  two  hours  of  CPU  time.  Parametric  studies  were  performed  to  gain 
fundamental  insight  into  process-induced  stress  development  in  thick-section  thermoset  laminates. 
Results  are  presented  and  discussed  in  the  next  section. 

3  Results  and  Discussion 

In  this  section  input  for  the  process-induced  stress  model  is  first  summarized.  Predictions  of 
the  effective  composite  in-plane  modulus  and  process-induced  strain  development,  associated  with 
the  cure  reaction  of  the  resin,  are  presented.  Several  aspects  of  the  analysis  are  then  validated 
by  comparison  with  experimental  results  of  modulus  and  laminate  curvature  found  in  previously 
published  literature.  Parametric  studies  are  presented  that  focus  on  the  effects  of  the  autoclave 
temperature  cure  cycle,  laminate  thickness,  stacking  sequence,  and  resin  cure  shrinkage  on  process- 
induced  stress  and  deformation  development. 
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Figure  3:  Process-induced  stress  modeling  flow  diagram. 
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3.1  Input  Parameter  Summary 


Required  model  input  includes  thermal  properties,  kinetic  parameters  and  mechanical  properties. 
This  input  for  glass/polyester  and  graphite/epoxy  was  summarized  in  Tables  1  through  4  in  the 
previous  section.  Typical  autoclave  temperature  cure  cycles  for  glass/polyester  and  graphite/epoxy 
employed  in  the  many  of  the  parametric  studies  presented  in  this  work  are  illustrated  in  Figure  4. 
Additional  input  (specified  where  appropriate)  includes;  total  laminate  thickness,  stacking  sequence 
and  volumetric  resin  shrinkage. 

As  with  all  numerical  solution  techniques,  the  potential  for  computational  round-off  error  exists. 
Convergence  studies  were  conducted  to  evaluate  the  sensitivity  of  the  analysis  to  nodal  spacing  and 
time  step  size  increment.  A  nodal  spacing  of  2.30xl0-3  m  and  a  time  step  of  1.0  second  were  found 
to  yield  converged  solutions  for  laminates  up  to  7.62  cm  thick  and  were  used  to  generate  all  results 
presented  in  this  investigation. 

3.2  Composite  Behavior  During  Cure 

3.2.1  Modulus 

The  influence  of  resin  modulus  on  the  effective  composite  properties  during  cure  is  quantified 
through  the  micromechanics  model.  Effective  transverse  and  in-plane  shear  moduli  of  glass/polyester 
are  presented  in  Figure  5  under  isothermal  cure  at  100°C.  The  rapid  increase  in  instantaneous 
moduli  is  indicative  of  the  rapid  curing.  Gradients  in  this  type  of  material  response  represent  a 
significant  mechanism  for  inducing  stress  in  thick-section  thermosetting  laminates  during  the  curing 
process. 


3.2.2  Process-Induced  Strains 

Volumetric  shrinkage  of  the  resin  is  occurring  simultaneously  with  effective  mechanical  property 
stiffening  effects.  Micromechanics  predictions  of  the  effective  macroscopic  composite  chemical 
shrinkage  strain  account  for  the  combined  influence  of  the  volumetric  resin  shrinkage  and  instan- 
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Figure  4:  Typical 


TIME  (MINUTES) 


ave  cure  cycles  for  glass/polyester  and  graphite/epoxy. 
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Modulus  (MPa) 


Figure  5:  Giass/polyester  moduli  during  isothermal  cure  at  100  °C. 
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taneous  resin  properties.  The  shrinkage  strain  in  glass/polyester  under  isothermal  cure  at  1Q0°C 
is  illustrated  in  Figure  6  for  a  total  volumetric  resin  shrinkage  of  vjh=6%.  While  the  longitudinal 
component  is  small,  due  to  the  constraining  fiber  stiffness,  the  effective  transverse  component  is 
significant,  exceeding  1%.  Gradients  in  these  shrinkage  strains  also  represent  a  significant  source 
for  stress  development  not  recognized  in  traditional  residual  stress  analyses. 

3.3  Comparison  With  Experimental  Results 

Hahn  and  Kim  subjected  several  graphite/epoxy  (T300/3501-6)  unidirectional  and  [0/90]  unsym- 
metric  laminates  to  interrupted  cure  cycles  to  monitor  the  development  of  transverse  modulus  (in 
the  unidirectional  laminates)  and  curvature  induced  by  cure  shrinkage  and  thermal  strains  (in  the 
unsymmetric  laminates)  [36].  Laminates  were  approximately  1  mm  in  thickness  and  the  cure  cy¬ 
cles  followed  the  manufacturer’s  recommended  cure  cycle  (see  Figure  4).  At  each  of  a  series  of 
points  during  cure,  a  unidirectional  and  a  cross-ply  laminate  were  cooled  to  room  temperature  at 
3°C/min.  The  unidirectional  laminates  were  tensile  tested  to  determine  the  transverse  modulus, 
and  the  deflections  of  the  unsymmetric  laminates  were  measured  to  determine  process-induced 
curvatures.  Figure  7  shows  excellent  agreement  between  the  measured  and  predicted  transverse 
modulus  development  based  on  the  material  model  presented  herein.  Note  that  the  horizontal  axis 
refers  to  time  of  cure  stop.  This  indicates  the  time  at  which  the  cure  cycle  was  interrupted  and 
cooling  to  room  temperature  initiated,  and  does  not  indicate  the  time  of  the  full  cycle  to  which  the 
laminate  was  subjected  as  not  all  laminates  were  cooled  from  the  same  temperature. 

Figure  8  shows  the  development  of  non-dimensionalized  curvature  (with  respect  to  laminate 
thickness)  as  a  function  of  the  time  of  cure  stop  for  total  epoxy  resin  volumetric  shrinkages  of  1.0 
and  1.5  %.  Again  we  see  good  agreement  between  measured  and  predicted  values.  Significant 
over-prediction  of  the  curvature  early  in  the  cycle  is  thought  to  arise  from  the  fact  that  the  highly 
viscoelastic  behavior  of  the  resin,  which  has  not  yet  reached  gelation,  is  not  incorporated  in  the 
present  model.  Also,  the  fracture  surfaces  of  the  tensile  tested  specimens  in  the  region  of  the 
curvature  over-prediction  exhibited  little  or  no  interfacial  bonding  between  the  fibers  and  the  resin, 
leading  to  a  relatively  unconstrained  resin  and  reduced  stress  development.  By  defining  the  gel 
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Figure  6:  Glass/polyester  chemical  shrinkage  strain  during  isothermal  cure  at  100 
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Figure  7:  Experimental  vs  predicted  composite  transverse  modulus  development  in  a  unidirectional 
graphite/epoxy  laminate. 
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Figure  8:  Experimental  vs  predicted  curvature  during  cure  of  an  unsymmetric  [90/0]  graphite/epoxy 
laminate. 


26 


point  as  the  point  at  which  the  viscosity  of  the  resin  reaches  100  Pa-Sec,  Hahn  and  Kim  note 
that  the  curvature  and  the  transverse  modulus  both  increase  sharply  immediately  after  the  gel 
point,  possibly  indicating  that  the  development  of  residual  stress  can  be  significantly  affected  by 
the  mechanical  properties  of  the  resin  before  full  cure  has  been  reached.  Good  agreement  with 
the  experimental  results  provides  confidence  in  the  accuracy  of  the  analysis  for  predicting  stress 
development  during  the  curing  process. 

3.4  Thickness  Effects 

To  gain  an  appreciation  for  the  mechanisms  introduced  in  this  investigation,  complex  temperature 
and  degree  of  cure  distributions  are  illustrated  in  glass/polyester  laminates  of  varying  thickness 
exposed  to  the  same  cure  cycle.  The  glass/polyester  cure  cycle,  indicated  in  Figure  4,  yields  a 
maximum  temperature  exotherm  in  a  2.54  cm  thick  laminate  at  164  minutes.  Figure  9  illustrates 
the  influence  of  laminate  thickness  on  temperature  distributions  at  164  minutes  into  this  cure  cycle. 
Temperature  exotherm  increases  with  laminate  thickness  below  2.54  cm  while  the  5  08  cm  laminate 
is  still  heating.  Corresponding  degree  of  cure  distributions  are  illustrated  in  Figure  10.  Degree 
of  cure  gradients  are  most  severe  at  this  point  in  the  2.54  cm  laminate.  The  interior  of  the  2.54 
cm  laminate  is  essentially  cured.  In  contrast,  the  5.08  cm  laminate  is  essentially  uncured  at  the 
interior.  At  a  later  point  in  the  cure  cycle,  the  thicker  laminate  will  exotherm  and  its  distributions 
will  reverse  shape  similar  to  the  thinner  laminates,  ultimately  developing  more  severe  gradients. 
These  results  demonstrate  the  complex  temperature  and  degree  of  cure  gradients  which  develop 
during  the  curing  process  of  thick  thermoset  laminates.  These  gradients  in  turn  have  a  profound 
influence  on  the  evolution  process-induced  stress  and  deformation. 

Residual  process-induced  transverse  stress  distributions  in  1.38,  1.85  and  2.54  cm  thick  unidirec¬ 
tional  glass/polyester  laminates  were  predicted  for  the  cure  cycle  shown  in  Figure  4  with  1^=6%. 
The  results  are  presented  in  Figure  11.  These  significant  stress  distributions  remain  after  complete 
cure  and  uniform  ambient  temperature  conditions  within  the  laminate  exist.  The  parabolic  stress 
profiles  are  self-equilibrating  tensile  stresses  at  exterior  and  compressive  stresses  at  the  interior 
regions.  This  profile  is  unique  to  laminates  that  cure  predominantly  from  the  inside  to  the  outside 
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Figure  9:  Temperature  distributions  in  glass/polyester  laminates  at  164  minutes. 
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Figure  10:  Degree  of  cure  distributions  in  glass/polyester  laminates  at  164  minutes. 
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DIMENSIONLESS  THICKNESS 


Figure  11:  Residual  process-induced  in-plane  transverse  stress  distributions  (£  <  2.54  cm). 
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surface  of  the  part.  Interior  regions  curing  first,  due  to  the  exotherm,  are  less  constrained  by  the 
uncured  exterior  region  during  shrinkage.  As  the  exterior  cures  and  contracts,  it  is  constrained  by 
the  cured  interior  region  of  appreciable  stiffness  and  the  process-induced  residual  stresses  shown 
in  Figure  11  develop.  Stress  gradients  are  accentuated  with  increasing  laminate  thickness.  For 
comparison  purposes,  it  is  pointed  out  that  stress  analyses  based  upon  the  traditional  sources  of 
residual  stress  associated  with  laminate  stacking  sequence  predict  that  unidirectional  laminates  are 
stress  free.  This  is  clearly  not  the  case  with  thick  laminates. 

Residual  process-induced  transverse  stress  distributions  in  unidirectional  glass/polyester  lami¬ 
nates  2.54,  5.08  and  7.62  cm  thick,  utilizing  the  same  cure  cycle  and  resin  shrinkage,  are  illustrated 
in  Figure  12.  A  dramatic  parabolic  stress  reversal  of  the  distributions  in  the  5.08  and  7.62  cm 
laminates  are  predicted  where  the  interior  regions  are  in  tension  and  the  exterior  regions  are  in 
compression.  This  reversal  results  from  the  change  to  an  outside  to  inside  cure  history  in  the  thicker 
laminates  and  clearly  demonstrates  the  significant  influence  the  laminate  thickness  has  on  the  evo¬ 
lution  of  process-induced  stress.  The  magnitude  of  peak  stresses  and  gradients  which  develop  are 
significant  enough  to  initiate  transverse  cracks  [3, 4, 5, 6]. 

3.5  Autoclave  Temperature  Cure  Cycle  Effects 

The  influence  of  the  autoclave  temperature  cure  cycle  history  on  stress  development  is  investi¬ 
gated  by  considering  the  curing  of  a  unidirectional  2.54  cm  glass/polyester  laminate  exposed  to 
the  temperature  histories  illustrated  in  Figure  13  and  a  resin  shrinkage  of  %.  Resulting 

residual  process-induced  transverse  stress  distributions  for  temperature  ramps  of  0.25,  0  5  and  1.0 
(°C/minute)  are  presented  in  Figure  14.  The  stress  distributions  developing  under  the  temperature 
ramps  of  0.25,  0-5  and  1.0  (°C/minute)  induce  compressive  stresses  at  the  interior,  while  the  more 
rapid  curing  process  results  in  the  development  of  tensile  interior  stresses.  These  results  demon¬ 
strate  the  importance  of  the  processing  history  on  residual  stress  development  in  thick  laminates. 
The  slower  thermal  ramps  create  an  inside  to  outside  cure  history  and  consequently  develop  inter¬ 
nal  compressive  stresses.  As  the  ramp  is  increased  the  cure  history  changes  to  an  outside  to  inside 
cure,  resulting  in  a  reverse  in  the  parabolic  stress  distribution. 
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DIMENSIONLESS  THICKNESS 


Figure  12: 


Residual  process-induced  in-plane  transverse  stress  distributions  (£  >  2.54  cm). 


32 


DIMENSIONLESS  THICKNESS 


Figure  14:  Glass/polyester  autoclave  temperature  cure  cycle  ramp  effects  on  transverse  stress 
distributions  ( l  =  2.54  cm). 
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3.6  Cure  Shrinkage  Effects 


The  influence  of  resin  shrinkage  on  residual  stress  development  was  also  investigated.  Total  vol¬ 
umetric  resin  shrinkages,  of  0,  1,  3  and  6  %  are  assumed  in  the  simulations  of  a  2.54  cm 

unidirectional  glass/polyester  laminate  cured  according  to  the  cure  cycle  in  Figure  4.  Residual 
transverse  stress  distributions  are  plotted  against  dimensionless  thickness  in  Figure  15. 

The  plots  indicate  that  the  magnitude  of  the  resin  shrinkage  significantly  influences  the  parabolic 
stress  distribution.  The  profile  for  vjh=0%  isolates  the  effect  of  non-uniform  curing  without  resin 
shrinkage.  It  is  interesting  to  note  that  for  vjh  =  l%  the  laminate  is  essentially  stress  free.  In  this 
case  the  processing  strains  due  to  chemical  shrinkage  are  effectively  equal  in  magnitude  but  opposite 
in  sign  to  the  thermal  expansion  processing  strains. 

3.7  Unsymmetric  Curing  Effects 

All  simulations  presented  to  this  point  exhibited  symmetric  curing  about  the  mid-plane  of  the 
laminate.  The  processing  histories  imposed  on  the  laminates  induced  either  inside-to-outside  or 
outside-to-inside  curing  patterns.  Consequently,  all  residual  stress  profiles  were  shown  to  be  sym¬ 
metric  about  the  laminate  mid-plane.  When  the  thermal  history  is  unsymmetric,  the  resulting 
stress  profiles  will  also  be  unsymmetric.  The  following  example  demonstrates  the  influence  of 
unsymmetric  curing  on  residual  stress  development. 

Unidirectional  glass/polyester  laminates  of  1.38,  2.54  and  3  39  cm  thickness  were  subjected  to 
unsymmetric  boundary  conditions.  The  autoclave  temperature  cure  cycle  illustrated  in  Figure  4  was 
prescribed  on  the  bottom  surface  of  the  laminate.  An  insulated  boundary  condition  was  imposed 
on  the  top  surface.  A  volumetric  resin  shrinkage  of  u^=6%  was  assumed.  In  these  simulations, 
the  cure  sweeps  from  the  bottom  to  the  top  surfaces.  The  resulting  residual  transverse  stress 
distributions  are  illustrated  in  Figure  16.  Similar  gradients  are  predicted  for  longitudinal  stress 
distributions.  Note  that  the  stress  profiles  are  self-equilibrating.  It  is  interesting  to  contrast  the 
magnitude  of  the  peak  transverse  stress  in  the  2.54  cm  laminate  presented  here  to  that  which 
developed  in  the  symmetrically  cured  laminate  presented  in  Figure  11.  The  symmetric  curing 


35 


process  resulted  in  transverse  stresses  ranging  between  -8.5  and  12.0  MPa  while  the  unsymmetric 
curing  process  reduced  stresses  substantially  (-2.0  and  2.0  MPa).  This  result  is  not  surprising  since 
the  unsymmetric  cure  induces  warpage  in  the  laminate,  reducing  the  magnitude  of  the  in-plane 
residual  stress. 

3.8  Stacking  Sequence  Effects 

The  parametric  studies  presented  to  this  point  have  focused  on  the  unidirectional  laminate  config¬ 
uration.  This  enabled  the  mechanisms  of  spatial  solidification  to  be  isolated,  without  the  influence 
of  the  traditionally  recognized  stacking  sequence  effects.  In  the  following  examples,  the  influence 
of  processing  on  stress  development  in  [0/90]  glass/polyester  cross-ply  laminates  is  quantified. 

In  the  first  example,  a  2.54  cm  laminate  of  [0/90]  symmetric  construction  with  8  repeating  units, 
([0/90]8»),  was  cured  symmetrically  according  to  the  glass/polyester  cure  cycle  with  v^=6.0%  (see 
Figure  4).  The  resulting  tranverse  stress  distributions  are  shown  in  Figure  17.  Such  stress  profiles 
are  similar  to  what  would  be  predicted  by  traditional  laminated  plate  theory  analysis  that  assumes 
a  stress-free  temperature.  In  this  simulation,  spatial  solidification  is  present  but  negligible  and  one 
concludes  that  stacking  sequence  effects  dominate. 

In  thicker  laminates,  however,  spatial  solidification  cannot  be  neglected.  Consider  the  resulting 
stress  profiles  in  symmetric  [0/90]„4  laminates  with  n=8,  16  and  24  (corresponding  to  total  laminate 
thicknesses  of  2.54,  5.08  and  7.62  cm,  respectively)  shown  in  Figure  18.  The  profiles  represent  the 
envelope  of  peak  normal  stresses  in  the  0°  or  90°  plies  through  the  laminate  thickness.  The  stresses 
are  discontinuous  across  ply  interfaces  and  are  equal  and  opposite  in  magnitude.  Processing  is 
shown  to  substantially  increase  the  stress  gradient  with  increasing  laminate  thickness.  In  contrast 
to  the  2.54  cm  laminate,  the  process-induced  stress  in  the  thicker  laminates  exhibit  significant 
gradients  through  the  thickness  that  are  also  much  greater  in  magnitude  in  comparison  with  the 
traditional  stacking  sequence  stress.  This  result  further  illustrates  the  significant  influence  the 
processing  history  has  on  stress  development  in  thick-section  thermosets,  even  for  laminates  of 
arbitrary  stacking  sequence. 
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Figure  16:  Influence  of  unsymmetric  curing  on  transverse  residual  stress  profiles  in  glass/polyester 
laminates. 
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Figure  17:  Transverse  residual  stresses  in  a  [0/90)8*  laminate. 
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Figure  18:  Normal  residual  stress  envelopes  in  [0/90]ni  laminates. 
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3.9  Stress-Free  Temperature 


Traditional  residual  stress  predictions  are  based  on  a  uniform  stress-free  temperature,  commonly 
assumed  to  be  the  cure  temperature.  The  uniform  temperature  drop  to  the  operating  temperature 
defines  the  temperature  difference  used  to  compute  residual  stresses.  In  the  following  example, 
the  influence  of  resin  shrinkage  on  the  assumed  stress-free  temperature  of  a  [Q/90]3  glass/polyester 
laminate  is  quantified. 

The  laminate  thickness  is  chosen  sufficiently  thin  to  minimize  temperature  and  degree  of  cure 
gradients  during  processing.  In  this  case,  stacking  sequence  is  the  dominant  source  of  residual 
stress.  The  autoclave  cure  cycle  is  defined  in  Figure  4.  Ply  stresses  are  predicted  for  volumetric 
resin  shrinkages  between  0  and  6%.  An  equivalent  stress-free  temperature  is  defined  from  our 
predictions  for  comparison  with  the  traditional  methodology  based  on  the  assumption  that  A T=- 
106°C  (the  difference  between  ambient  and  the  cure  temperature).  A  dimensionless  stress-free 
temperature  difference  is  defined  as  the  ratio  of  the  predicted  to  the  assumed  temperature  change. 
In  Figure  19,  dimensionless  stress-free  temperature  is  plotted  against  the  volumetric  resin  shrinkage. 
Only  for  a  volumetric  shrinkage  of  1%  is  the  traditional  assumption  valid  (AT =-106°C).  The  results 
demonstrate  the  significant  influence  that  resin  shrinkage  can  have  on  the  assumed  stress- free 
temperature,  and  thus  the  magnitude  of  the  resulting  residual  stress  distributions. 

Another  point  worth  noting  is  that  in  our  simulation,  the  evolution  of  stress  throughout  the 
curing  process  is  predicted  and  therefore  no  stress-free  temperature  assumption  is  required.  Fur¬ 
thermore,  the  severe  gradients  in  temperature  and  degree  of  cure  that  were  shown  to  develop  in 
thick  laminates  suggest  that  the  traditional  approach  is  not  appropriate.  This  is  apparent  from 
the  fact  that  significant  residual  stresses  were  shown  to  develop  in  thick  unidirectional  laminates. 
Recall  that  the  traditional  analysis  predicts  that  unidirectional  laminates  are  stress  free.  Therefore, 
the  traditional  residual  stress  analyses  are  not  appropriate  for  thick-section  laminates  in  general. 
The  appropriateness  of  the  assumption  depends  on  many  factors  such  as  laminate  thickness,  resin 
shrinkage  and  the  processing  history. 
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Dimensionless  Stress  Free  Temperature 


Volumetric  Resin  Shrinkage  (%) 


Figure  19:  Influence  of  resin  shrinkage  on  dimensionless  stress-free  temperature. 
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4  Conclusions 


A  fundamental  study  of  process-induced  stress  and  deformation  in  thick-section  thermosetting 
composite  laminates  was  presented.  A  methodology  was  proposed  for  predicting  residual  stress  de¬ 
velopment  during  the  curing  process  that  does  not  require  a  stress-free  temperature  to  be  assumed. 
A  one-dimensional  cure  simulation  analysis  was  coupled  to  an  incremental  laminated  plate  theory 
model  to  study  the  relationships  between  complex  gradients  in  temperature  and  degree  of  cure, 
and  process-induced  residual  stress  and  deformation. 

Material  models  were  proposed  to  describe  the  mechanical  properties,  thermal  and  chemical 
strains  of  the  thermoset  resin  during  cure.  The  material  models  were  incorporated  into  a  microme¬ 
chanics  model  to  predict  the  effective  mechanical  properties  and  process-induced  strains  of  the 
composite  during  cure.  Thermal  expansion  and  cure  shrinkage  contribute  to  changes  in  material 
specific  volume  and  were  shown  to  represent  important  sources  of  internal  loading  included  in  the 
analysis. 

Temperature  and  degree  of  cure  gradients  that  develop  during  the  curing  process  induce  me¬ 
chanical  property  and  process-induced  strain  gradients  in  the  composite  laminate.  These  gradients 
represent  important  mechanisms  that  contribute  to  stress  and  deformation  development  not  previ¬ 
ously  considered  in  traditional  residual  stress  analyses  of  laminated  composites.  Mode!  predictions 
of  cure  dependent  epoxy  modulus  and  curvature  in  unsymmetric  graphite/epoxy  laminates  were 
found  to  be  in  good  agreement  with  experimental  data  presented  in  previously  published  literature. 

The  effects  of  processing  history  (autoclave  temperature  cure  cycle),  laminate  thickness,  resin 
cure  shrinkage  and  laminate  stacking  sequence  on  the  evolution  of  process-induced  stress  and  defor¬ 
mation  in  thick-section  glass/polyester  and  graphite/epoxy  laminates  during  cure  were  quantified. 
Processing  of  unidirectional  laminates  was  investigated  to  provide  fundamental  insight  into  the  im¬ 
portant  mechanisms  governing  residual  stress  development  during  cure.  While  traditional  sources 
of  residual  stress  associated  with  laminate  stacking  sequence  predict  that  unidirectional  laminates 
are  stress  free,  our  results  clearly  indicate  the  potential  for  significant  process-induced  stress  devel¬ 
opment  in  thick  thermoset  laminates. 
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Residual  stresses  were  shown,  in  general,  to  increase  with  increasing  laminate  thickness,  au¬ 
toclave  temperature  ramp  and  resin  shrinkage.  Unsymmetric  curing  was  shown  to  significantly 
reduce  residual  stress  in  comparison  with  a  symmetric  curing  process.  Spatial  solidification  and 
stacking  sequence  effects  were  found  to  be  equally  important  in  thick-section  cross-ply  laminates. 
Furthermore,  the  results  indicated  that  the  assumption  of  a  stress-free  temperature  for  residual 
stress  predictions  is  not  appropriate  for  thick-sections  in  general. 

The  present  study  has  demonstrated  fundamental  mechanisms  that  contribute  to  residual  stress 
development  that  have  not  been  previously  considered  in  traditional  residual  stress  analyses  of  lam¬ 
inated  composites.  The  significant  magnitude  of  process-induced  residual  stress  identified  in  this 
work  substantiates  the  prevalent  concerns  associated  with  processing  difficulties  of  thick  thermoset¬ 
ting  laminates,  such  as  delaminations  and  matrix  cracking.  The  results  clearly  indicate  that  the 
mechanics  and  performance  of  thick-section  thermoset  laminates  are  strongly  dependent  on  pro¬ 
cessing  history. 
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A  Continuous  Fiber  Micromechanics  Model 


The  doubly  embedded  self-consistent  field  model  approach,  based  on  the  work  of  Hill  [37],  offers  a 
more  realistic  approximation  of  the  effective  lamina  properties  than  mechanics  of  materials  models 
based  on  the  uniform  stress  or  strain  assumption.  The  model  also  offers  mathematically  tractable 
closed-form  expressions  for  the  effective  lamina  mechanical  properties  and  expansion  strains.  The¬ 
oretical  background  of  the  model  is  presented  in  greater  detail  elsewhere  [32]. 


In  the  equations  presented  below,  the  subscripts  1,  2,  and  3  refer  to  the  directions  in  the  principle 
coordinate  system  of  the  lamina.  Subscripts  m  and  /  correspond  to  the  matrix  (resin)  and  fiber 
properties,  respectively.  Fiber  volume  fraction  of  the  lamina  is  denoted  as  vj  and  k  is  the  isotropic 
plane  strain  bulk  modulus  defined  by 


k  = 


E 

2(1  -v-2v*) 


(24) 


A.l  Engineering  Constants 


The  following  equations  define  the  transversely  isotropic  engineering  constants  of  the  lamina.  The 
longitudinal  Young’s  modulus: 

[4(^12 m  ~  ^/)fc/fcmG23m(l  ~  Vf)Vf 


E\  =  Eijvj  4-  f?im(l  -  Vj)  + 


[(£/  -t-  G23m)km  +  (kj  -  km)Gnmvf\ 


The  major  Poisson’s  ratio: 

*>12  =  ^13  =  VUfvJ  +  y12m(l  -  Vf)  + 
The  in-plane  shear  modulus: 

Gl2  —  G 13  =  Gl2m 

The  transverse  shear  modulus: 


(^12m  -  Vl2f)(km  -  fc/)C?23m(l  ~  vf)vJ 

(kf  +  C?23m  )km  +  (kj  —  km  )(?23m  Vy 


(Gl2f  +  Gl2m)  +  (Gi2f  ~  G\2m)vf 
(Guf  +  Gi2m)  -  {G\2f  ~  Gi2m)f/ 


q  _  C23m[fcm(G23m  +  *-*23/)  ~t~  2G23/(?23m  +  km(G23f  ~  C23m)n/] 
km(G23m  +  G23/)  +  2G23/C23m  ~  {km  -(-  2G23m)(G23/  ~  G23m)Vf 


(25) 


(26) 


(27) 


(28) 
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The  transverse  Young’s  modulus: 


E 2  E3  (l/4tj)  +  (I/4G23)  +  ivn/ E\) 

where  k?  is  the  effective  plane  strain  bulk  modulus  of  the  composite  given  by 

_  (*/  +  G23m)km  +  ( kf  -  km)G23mVf 
T  {kf  +  G23m)  -  (*7  -  km)vj 


The  transverse  Poisson’s  ratio: 


^23  = 


2Eikj  —  £i£2  —  E2 

2E\kj 


(29) 


(30) 


(31) 


A. 2  Expansional  Strains 


Expansional  strain  expressions  are  used  to  evaluate  the  in-plane  chemical  shrinkage  strains  and 
the  effective  coefficients  of  thermal  expansion  of  the  lamina  (aj  =  and  =  *2)-  The  following 
equations  define  the  transversely  isotropic  expansional  strains  of  the  lamina.  The  longitudinal 
direction  expansional  strain: 

(\fE\fv]  +  ZlmE\m{\  ~  Vf) 


tl  = 


E\fVf  +  Elm(l  -  Vf) 
The  transverse  direction  expansional  strain: 


(32) 


£2  =  f3  =  («2/  +  "12/fl/W  +  (£2m  +  ^12mflm)(l  ~  Vf)  - 

wnm + -w>  -  v/)i  [■i'^:r.Ei',-i:>r,)] 


(33) 
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